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ABSTRACT 


This  dissertation  investigates  the  plane  elasticity  problem 
of  a finite,  rigid  rectangular  block  partially  embedded  in,  and 
perfectly  bonded  to  an  elastic  half  space.  The  bond  thickness  is 
assumed  to  be  sufficiently  thin,  so  that  there  is  no  discontinuity 
in  the  displacements  of  the  bonded  surfaces. 

The  problem  is  formulated  by  superposition  of  the  solutions 
to  the  problems  of  horizontal  and  vertical  line  inclusions  beneath 
an  elastic  half  space,  which  are  derived  from  integral  transform 
techniques.  Substitution  of  these  results  into  the  boundary  condi- 
tions appropriate  for  the  embedded  block  problem  leads  to  a system 
of  six  coupled  singular  integral  equations,  whose  unknowns  are  the 
normal  and  shear  stress  discontinuities  between  the  bonded  surfaces. 

Three  distinct  sets  of  loads  are  applied  to  the  embedded 
block,  so  that  it  either  translates  without  rotation  in  the  x-  or 
y-  direction,  or  rotates  about  an  axis  in  the  z-direction.  In  all 
three  cases,  by  taking  advantage  of  the  symmetric  nature  of  the 
problems,  the  number  of  the  governing  singular  integral  equations 
is  reduced  from  six  to  four.  These  equations  may  be  solved  numeri- 
cally by  employing  the  collocation  scheme  introduced  by  Erdogan, 
Gupta,  and  Cook,  provided  that  they  be  supplemented  with  four  sub- 
sidiary conditions.  The  requirement  for  additional  conditions  is 
dictated  by  the  order  of  the  singularities  at  the  corners  of  the 
block,  established  through  an  asymptotic  expansion  of  the  singular 
integral  equations  in  the  vicinity  of  these  corners. 


iv 


The  problem  of  two  vertical  inclusions  perpendicular  to  the 
free  surface  of  the  half  space,  and  loaded  in  a manner  analogous 
to  the  embedded  block  problem,  is  also  formulated  for  each  case. 
These  problems,  besides  their  also  being  of  interest,  provide  a fur- 
ther estimate  of  the  accuracy  of  the  numerical  results. 

The  numerical  analysis  was  carried  out  for  various  geometries 
for  the  case  of  plane  strain.  Several  important  physical  quantities 
are  computed,  such  as  the  diffusion  of  the  load  from  the  block  into 
the  elastic  half  space,  and  the  local  stress  distribution  around 
the  block,  excluding  the  singular  points. 
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CHAPTER  I 


INTRODUCTION 

This  dissertation  is  concerned  with  the  stress  analysis  of  an 
elastic  half  space,  in  which  a perfectly  bonded,  rigid  rectangular 
block  is  partially  embedded.  The  bond  thickness  is  assumed  to  be 
sufficiently  thin  so  that  the  displacements  of  the  bonded  surfaces 
are  continuous.  Within  this  context  the  problem  is  one  of  plane 
strain  and  can  be  considered  appropriate  to  calculations  that  involve 
the  stress  distributions  around  foundations  where  the  out-of-plane 
dimensions  are  very  large  when  compared  with  the  length  or  width  of 
the  rectangle.  Alternatively,  the  problem  can  be  viewed  as  one  of 
generalized  plane  stress  analyzing  the  load  diffusion  from  a finite, 
rigid  rectangular  insert,  partially  embedded  within  a semi-infinite 
sheet,  where  the  axis  of  the  insert  is  perpendicular  to  the  edge  of 
the  sheet.  The  geometry  and  coordinate  system  for  such  a block  is 
shown  in  Fig.  1,  where  y is  the  shear  modulus  and  K is  related  to 
Poisson's  ratio  by  K * 3 - 4v  (plane  strain)  or  K > (3-v)/(l+v)  (plane 
stress),  and  v is  Poisson's  ratio. 

Problems  such  as  the  one  considered  here  have  had  a long 

history,  which  is  best  summarized  in  a paper  by  Muki  and  Sternberg 

[1] , who  consider  the  diffusion  of  load  from  a transverse  tension  bar 

into  a semi-infinite  elastic  sheet.  They  reconsider  the  problem 

initially  posed  by  Reissner  [2] , who  considered  the  load  transfer 

from  a transverse  stringer,  a finite  segment  of  which  overlaps  with, 
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and  is  continuously  bonded  to,  a semi-infinite  elastic  sheet.  Theirs 
is  essentially  a contact  problem  in  which  the  finite  stringer  is  at- 
tached to  the  semi-infinite  sheet  and  their  objective  is  to  obtain  a 
systematic  reduction  of  the  problem  to  a Fredholm  integral  equation 
via  two  contact  models:  i.e.  line-contact  or  area-contact. 

The  present  analysis  differs  from  their*  in  that  the  sheet  is 
assumed  to  be  cut  so  that  it  has  a finite  rectangular  notch  at  its  sur- 
face. A perfectly  matching  rigid,  rectangular  insert  is  then  bonded 
within  this  notch. 

Loading  is  applied  to  this  embedded  insert  so  that  it  either 
translates  without  rotation  in  the  vertical  or  horizontal  direction, 
or  rotates  due  to  an  applied  moment.  The  plane  strain  case  corres- 
ponds to  an  infinite,  rectangular  block  embedded  within  an  elastic 
half  space  undergoing  vertical  or  horizontal  translation,  or  rotation. 

It  should  be  noted  that  Muki  and  Sternberg  have  also  dealt  with 
the  three-dimensional  problem  of  load  diffusion  from  an  axially 
loaded  rod  to  a half  space  [3,4],  In  [ 3 J the  problem  of  the  axial 
force  decay  in  an  infinite  cylindrical  elastic  bar  bonded  bo  an  infinite 
medium  is  dealt  with.  An  approximate  solution  scheme  for  cross- 
sections  of  arbitrary  shape  is  developed  and  compared  to  the  exact 
one  for  the  case  of  a circular  bar.  This  scheme  is  then  used  in  [4] 
to  solve  the  problem  of  load  diffusion  from  a bar  of  arbitrary  uni- 
form cross-section  partially  embedded  in,  and  bonded  to  a semi- 
infinite solid. 

The  solution  of  the  considered  problem  is  developed  in  succes- 
sive stages.  In  Chapter  II  Green's  functions  are  established  for  a 
horizontal  and  vertical  line  load  beneath  an  elastic  half  space. 
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This  formulation  allows  the  boundary  value  problem  to  be  written  as 
having  boundary  conditions  only  on  the  surface  of  the  notch.  The 
results  are  derived  using  integral  transform  techniques,  but  they  can 
also  be  derived  by  other  methods.  The  primary  references  used  for 
these  derivations  were  the  book  by  Sneddon  [5] , and  the  tables  edited 
by  Erdelyi  [6]. 

Chapter  III  deals  with  the  vertical  translation  of  the  finite 
rigid  rectangular  block  partially  embedded  in,  and  perfectly  bonded 
to  an  elastic  half  space,  under  an  applied  load  2P  acting  in  the 
negative  y-direction.  The  appropriate  boundary  conditions  are  written 
by  superposing  the  horizontal  and  vertical  inclusion  results.  The 
boundary  conditions  become  the  data  for  a system  of  singular  integral 
equations,  which  can  be  solved  numerically.  It  is  noted  that  there 
are  four  points  on  the  block  that  are  singular:  the  intersection  of 
the  block  with  the  free  surface,  (+c,0),  and  the  lower  corners  of  the 
block,  (+c,h).  These  singularities  can  be  accounted  for  in  the  ana- 
lysis by  the  establishment  of  appropriate  corner  conditions  through 
an  asymptotic  expansion  of  the  governing  singular  integral  equations 
in  the  vicinity  of  these  singularities.  It  should  be  noted  here 
that  the  solution  of  the  transcedental  equation  derived  for  the  corners 
(+c,h)  leads  to  two  negative  roots.  Therefore,  they  both  lead  to 
singular  contributions.  The  singular  stress  field  is  dominated  by  the 
largest  singularity,  and  this  is  the  one  used  in  the  present  analysis. 
To  the  author's  knowledge,  the  only  work  that  considers  both  negative 
roots  at  a corner  such  as  the  one  encountered  here,  is  a paper  by 
Westmann,  [7],  which  deals  with  a wedge  bonded  to  a half  plane  along 
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a finite  length.  The  solution  for  that  problem  was  developed  using 
Mellin  transforms,  and  numerical  results  were  obtained  by  use  of  a 
modified  finite  difference  method.  However,  accounting  for  both  nega- 
tive roots  is  beyond  the  scope  of  this  dissertation.  Finally,  in 
order  to  obtain  a further  estimate  of  the  accuracy  of  the  numerical 
results,  an  analogous  problem  is  formulated.  This  is  the  problem  of 
two  rigid  vertical  inclusions  of  length  equal  to  the  sides  of  the 
rigid  block,  and  separated  by  a distance  equal  to  the  base  of  the 
block,  partially  embedded  in  and  bonded  to  an  elastic  half  space, 
each  loaded  by  P in  the  negative  y-direction. 

The  case  of  horizontal  displacement  of  the  rigid  block  loaded 
by  2Q,  acting  in  the  negative  x-direction,  and  a moment  M is  con- 
sidered in  Chapter  IV.  The  analogous  case  for  the  two  inclusion  prob- 
lem is  also  examined  in  the  same  chapter. 

Chapter  V investigates  the  rotation  of  the  rigid  block  under  the 
action  of  a moment  M.  The  two  inclusion  problem  subjected  to  analogous 
loading  is  also  formulated  in  that  chapter. 

In  Chapter  VI  the  method  for  obtaining  numerical  results  is 
discussed.  The  collocation  scheme  introduced  by  Erdogan,  Gupta,  and 
Cook  [8]  is  used  to  solve  the  system  of  singular  integral  equations 
governing  each  of  the  problems  formulated  in  Chapters  III-V.  The 
method  of  solution  allows  the  calculation  of  certain  inportant  physi- 
cal quantities  such  as  the  diffusion  of  load  from  the  block  into  the 
elastic  half  space,  and  also  the  local  stress  distribution  around  the 
block,  excluding  the  singular  points. 


CHAPTER  II 
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FORMULATION  AND  BASIC  EQUATIONS 


HORIZONTAL  INCLUSION  IN  A HALF  SPACE 

The  equations  to  be  obtained  here  are  for  a bonded,  rigid  inclu- 
sion of  length  2c  located  at  a depth  h beneath  the  free  surface,  y = 0, 
of  an  elastic  half  space.  The  inclusion  is  assumed  to  have  zero  thick- 
ness and  its  geometry  and  coordinate  system  are  shown  in  Fig.  2.  If 
the  discontinuities  in  the  normal  and  shear  stresses  are  designated  by 
A (x)  and  B (x) , respectively,  then  the  following  conditions  must  be 
satisfied: 


(2) 

- T(1) 

= A (x) 

yy 

yy 

(2) 

(i) 

=*  B (x) 

- T 

xy 

xy 

(2.1) 


In  addition  to  the  above  equations,  the  following  continuity 
and  boundary  conditions  must  be  satisfied: 


0 


. (2) 


.(1) 


y = h,  -c  < x < c 


(2.2) 


y y 

t (x,0)  * t (x,0)  “0,  0 <_  |x|  < “ (2.3) 

yy  xy 

where  the  superscripts  (1),  (2)  define  the  regions  above  and  below  the 
inclusion. 

In  terms  of  the  stress  discontinuities,  A(x)  and  B(x)  defined 
above,  the  displacement  derivatives  and  stresses  throughout  the  half 
space  are  given  as 
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The  functions  K ^ and  Kg^,  i = 1,2,. ..,7,  are  given  in  Appendix  A and 
they  are  rational  functions  of  x and  y with  both  numerators  and  denomina- 
tors being  polynomials  in  x and  y. 

Equations  (2.4)  - (2.10)  were  derived  using  integral  transform 
techniques  [5,6] , but  they  can  also  be  derived  by  other  methods. 


VERTICAL  INCLUSION  IN  A HALF  SPACE 


Here,  the  case  of  a vertical  inclusion  of  length  (b  - a)  perpen- 
dicular to  the  free  surface  of  an  elastic  half  space,  y * 0,  is  solved 
(Fig.  2) . Let  the  jumps  in  the  normal  and  shear  stresses  across  the 
inclusion  be  C(y)  and  D(y),  respectively;  then,  the  following  relation- 
ships must  holds 


,<*>  - T<‘>  - C(y) 


XX 


XX 


t(2)  - T(1)  - D(y) 
xy  xy 

U<2)  -U(1)  - 0 

X X 

„<*>  - „<»  - o J 

y y 


x ■ c,  a ^ y ^ b ; 


x * c,  a^y£b; 


T (X,0)  “ T (X,0) 

yy  *y 


0,  y - 0,  0 < |x|  < 


(2.11) 


(2.12) 

(2.13) 


For  this  case  the  superscripts  (1) , (2)  refer  to  the  regions  to  the 
left  and  right  of  the  inclusion,  respectively. 


In  terms  of  the  stress  jumps,  C(y)  and  D(y),  defined  by  Eqs. 
(2.11),  the  displacement  derivatives  and  stresses  within  the  half 
space  are  given  as 


. II 
3u 

,1} 

b 

,y ; t)dt} 

X 

3x 

(x,y) 

1 

4ir  (ic  + l)y 

U 

(x, 

,y;t)dt 

+ / 

D<t>ln 

(x, 

a 

a 

(2.14) 
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1 

4it(k  +1)u 

{/ 
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(x, 

a 

a 

(2.15) 
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3u 

b 
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,y;t)dt} 

X 
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X 
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a 
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W 

a 

D(t)Ss 

(x, 
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II 


du 


3y 


(*,y)  ■ 


4it(k  +1)u 


II , . 

Txx(x'y)  * 


II,  n 

V*>y) 


r+ty  4 c<t)Ins(x,yit)dt  + [ ■>(OLS5(x.yitwtl) 


2tt(k  +1) 


2tt  (*c  +1) 


TTiT  4 c<t)h„U.yit)dt  + J D<t)Ls6(x,y;t)dt)^ 


Ulhl)  4 c<t)LH7(x,y;t>dt  + { 


where  the  functions  and  Lg^,  i “ 1,2,..., 7,  are  given  in  Appendix  A. 
It  is  noted  that  the  superscript  II  is  used  to  identify  displacement 
derivatives  and  stresses  associated  with  the  vertical  inclusion  at  x = c, 
while  I was  used  in  connection  with  the  horizontal  inclusion  located 
at  y = h.  To  obtain  the  equations  associated  with  a vertical  inclusion 
of  length  (b  - a) , located  at  x = -c,  one  replaces  c by  -c  in  Eqs. 

(A. 15) - (A. 28) , C(t)  by  E(t),  and  D(t)  by  F(t)  in  Eqs.  (2 . 14) - (2 . 20) . 

For  this  case  the  superscript  III  will  be  used.  These  equations  will 
not  be  listed  here  but  will  be  incorporated  in  the  following  chapters. 


CHAPTER  III 


STRESS  ANALYSIS  OF  EMBEDDED  BLOCK:  VERTICAL  DISPLACEMENT 


VERTICAL  TRANSLATION  OF  RIGID  BLOCK 


The  results  of  the  preceding  chapter  are  now  superposed  to  formu- 
late the  problem  of  a rigid  block  perfectly  bonded  to  an  elastic  half 
space  and  loaded  as  shown  in  Fig.  1.  The  limits  of  integration  for  the 
vertical  inclusions  (solutions  II  and  III)  are  taken  as  a = 0 and 
b = h (Fig.  2).  The  boundary  conditions  for  this  problem  are  given 
next: 


du  3U1  3U11  3U111 

a -f  " if  (x’h)  + (x*h)  + IT-  (x*h)  " 0 


3u  du 

2.  . 2 

3x  3x 


3u  3u 


du 


II 


3u 


III 


(x,h)  + (x,h)  + — (x,h)  - 0 


, -c  < x < c,  y-h; 

(3.1) 


3u 


II 


3u 


III 


3y 

III 


(-c,y)  - 0 


3T"  3 f + 3 T"  ('C'y)  + 

3u  3u^  3u^  3u 

(_C’y)  + (_C»y)  + "3 (_C‘y)  “ ° 

3u  3u*  3u**  Su1** 

TT  m TT  (c,y)  + ~£~  (c.y)  + 


f,  0 < y < h,  x - -c  , 
(3.2) 


3y 

II 


3y 

III 


(c,y)  • 0 


3u  3u"  Su*"*-  3u 

(c*y)  + ij”  (c*y)  + (c’y)  " ° 


r*  0 < y <.  h»  x " c- 

(3.3) 


Equations  (3.1)- (3. 3)  represent  a system  of  six  singular  integral 
equations,  which  are  given  next: 
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c c h 

J Ats)!^  (x,h;s)ds  + / B(s)KSl(x,h;s)ds  + / C (t)LN1 (x,h;t)dt 

•c  -c  0 


h 

+ f D(t)L  (x,h;t)dt  + / E(t)iL  (x,h;t)dt  + / F(t)L  (x,h;t)dt  - 0, 
o o 1 0 

(3.4) 


h 

/ 

0 


h 

/ 

0 


-c  <_  x <_  c , 


v_  C •• 

/ A(8)KN2(x,h;s)ds  + / B(a)KS2(x,h;s)d8  + / C(t)LN2(x,h;t)dc 

_c 


+ j D(t)LS2(x,h;t)dt  + j E(t)L^2 (x,h;t)  + / F(t)LS2(x,h;t)dt  - 0 

-C  < X < C| 

c c h 

J A(a)KN3(-c,y;8)d8  + / B(s)KS3 (-c ,y;s)ds  + / CCOL^-c.y.tJdt 


(3.5) 


-c 


h h 

+ / D(t)Lc  (-c,y;t)dt  + / E(t)iL  (-c,y;t)dt  - 0,  0 < y < h>  (3.6) 

o »3  o 3 

c c h 

J A(8)KN4(-c,y;a)d8  + / B(s)KS4(-c,y;s)ds  + / C(t)LN4(-c,y;t)dt 

-c  -c  O' 

h h 

+ / D(t)L_  (-c,y;t)dt  + / F(t)Lg  (-c,y;t)dt  - 0,  0 < y < h,  (3.7) 

0 5,4  0 


c 

/ 

-C 


c c ** 

/ A^K^c.yjsJda  + / B(a)KS3 (c ,y;a)da  + / CCt)!^ (c ,y ;t)dt 


+ / E(t)LN3(c,y;t)dt  + / F(t)LS3(c,y;t)dt  - 0,  0 < y < h,  (3.8) 


C u ** 

J A(8)KN4(c,y;8)d8  + / B(s)KS4(c,y;s)ds  + / D(t)LS4(c,y;t)dt 


h h 

+ J E(t)LN4(c,y;t)dt  + / F(t)LS4(c,y;t)dt  - 0,  9 < y < h.  (3.9) 


The  barred  kernels  and  Lg^,  i ■ 1,2,..., 4,  are  readily  obtained 
from  the  corresponding  unbarred  kernels  given  in  Appendix  A by  replacing 
c by  -c.  It  is  also  noted  that  the  boundary  conditions  for  the  integral 
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equations  are  homogeneous.  However,  the  input  data  will  be  given  in 
the  form  of  subsidiary  conditions  that  involve  the  load  2P  applied  to 
the  block  and  the  constraints,  when  the  symmetry  appropriate  to  the  con- 
sidered problem  is  used. 

Since  Eqs . (3.4)- (3.9)  deal  with  the  geometry  of  corners,  their 
limit  near  the  corners  must  produce  relationships  which  hold  true  in 
the  neighborhood  of  such  corners.  These  relationships  may  be  directly 
obtained  from  the  governing  singular  integral  equations,  after  the 
contributing  parts  are  isolated  from  the  entire  equations.  First,  the 
order  of  the  singularity  at  the  intersection  of  the  free  surface  with  a 
vertical  inclusion  will  be  established  by  considering  the  appropriate 
contributory  terms.  Adjusting  the  coordinate  system  to  account  for  the 
right-hand  inclusion,  the  terms  leading  to  the  determination  of  the 
singularity  at  the  free  surface  are  given  as 


(3.10) 


h 

, K-+1  <i 

0 y-t  y+t  "(y+t) 

where  y is  approaching  zero.  A second  equation,  identical  to  (3.10),  is 
obtained  for  D(t).  Assuming  solutions  of  the  form, 


C (t)  = C(t)tn_1,  D (t)  * D(t)tn_1 


(3.11),  (3.12) 


where  0 < n < 1,  the  following  eigenvalue  equation  is  deduced  for  the 
determination  of  n: 


. 2 .irn,  (k+1)  2 1 2 

sin  l-f)  » — -n. 


(3.13) 


Except  for  a minor  change  in  notation,  Eq.  (3.13)  is  identical  to  that 
determined  by  Williams  [9]  for  a right-angle  corner  of  an  elastic  plate 
in  extension  with  fixed-free  boundary  conditions.  When  Eq.  (3.13)  is 
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solved  for  v » 0.3  (plane  strain),  the  result  is 

n * 0.71117.  (3.14) 

Conditions  at  the  corners  (c,h)  and  (-c,h)  can  be  established 
in  a similar  manner  through  the  use  of  Eqs.  (3.4)— (3.9) , where  the 
following  asymptotic  forms  are  assumed  for  the  sought  functions: 

A(s)  - A(s)(c2  -s2);  1 , B(s)  - B(s)  (c2  - s2)** 

C (t)  - C(t)tn_1  (h-t)C_1,...,F(t)  - F(t)tn"1(h-t);"1  (3.15) 

where  0 < 5 < 1.  The  corner  conditions  for  (c,h)  are  found  to  be 

k cos(ttc)B0  + (ic  - C)cos(^-)C0  + 5 sin(^-)D0  - 0 (3.16a) 

K cos(it5)A0  + c sin(^)C0  + (k  + ;)cos(^-)D0  - 0 (3.16b) 

S sin(^-)AQ  + (ic  + ;)cos(^-)B0  + k cos(irc)C0  - 0 (3.16c) 

(k  - Ocos(^)A0  + t sin(^-)B0  + k cos(its)D0  - 0,  (3.16d) 

and  for  (-c,h) 

ic  cos (tt^)BI  + (ic  - c)cos(^-)E0  - ? sin(-^-)F0  - 0 (3.17a) 

k cos(tt;)Aq  - i sin(^-)EQ  + (k  + ;)cos(^-)F0  - 0 (3.17b) 

(k  - C)cos(7t?)Aq  - 5 s in(^)Bg  + k cos(?r;)F0  - 0 (3.17c) 

sin(^-)Ag  + (k  + C)cos(^-)Bq  + k cos(irC)E0  - 0.  (3.17d) 

In  Eqs.  (3.16)  and  (3.17)  the  following  notation  is  used: 

A0  - A(c)(2c);'1/hn"1,  B0  - B(c)(2c)5‘1/hn‘1, 

A;  - A(-c)(2c)C“1/hn"1,  B!  - B(-c)(2c)5“1/hn"1, 

(3.18) 

C0»  D0«  E0*  F0  “ C(h),  D(h),  1(h),  F(h). 
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Equations  (3.16)  and  (3.17)  lead  to  the  same  eigenvalue  equation, 
namely 


[tc2sin2(^)  - £2]  [ic2sin2  - £2]  - 0, 

where  £ is  determined  from 

8in(^)  ■ 1 ~ c» 


(3.19) 


(3.20) 


which  is  in  agreement  with  Williams'  result  for  a fixed-fixed  3it/2  - 
corner  of  an  elastic  plate  in  extension  [9],  If  the  positive  sign 
in  Eq.  (3.20)  is  chosen  and  v = 0.3  (plane  strain),  the  result  is 

C = 0.59516  . (3.21) 

It  should  be  noted  that  Eq.  (3.20)  with  the  negative  sign  also  yields 
a root  in  the  interval  of  interest,  £ = 0.75904;  however,  since  the 
numerical  scheme  used  to  obtain  numerical  results  [8]  can  only  accom- 
modate one  of  these  roots,  the  gravest  one,  given  by  Eq.  (3.21)  will 
be  used. 

Equations  (3.4)- (3.9)  represent  a system  of  six  singular  integral 
equations,  which  along  with  appropriate  subsidiary  conditions  can  be 
solved  for  the  six  unknown  stress  discontinuities  A(s),  B(s),  ...»  F(t). 
However,  by  taking  advantage  of  the  symmetric  nature  of  the  problem, 
the  number  of  unknowns  can  be  reduced  to  four,  since  the  following 
relationships  must  hold: 

E(y)  = -C(y),  F (y)  = D(y).  (3.22) 

Thus,  the  subsidiary  conditions  associated  with  equations  (3.4)  — (3.9) 
are  not  needed, and  will  not  be  listed  here.  Equation  (3.22)  is  in- 
corporated to  write  Eqs.  (3.4)— (3.5) , and  (3.8)-(3.9)  in  the  following 
form: 
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c c 

/ A(s)KNi(x,h;s)ds  + / B/s) Kgl (x,h; s) ds 


h _ 

+ / C(t) [Ln1 (x,h;t)  - ln1 (x,h;t) ]dt 


h _ 

+ / D(t) (L  (x,h;t)  + L (x,hjt) ]dt  = 0,  -c  x < c ; (3.23) 

£ SI  SI 


c 

J A ( s ) K ^ ( x , h ; s ) ds 

-c 


c 

+ / B (s) (x,h; s) ds 
-c 


h _ 

+ / C (t) [Ln2 (x,h; t)  - Ln2  (x,h;t)]dt 


h 

+ / D (t)  [L  (x,h;  t)  + L_,  (x,h;t)  ]dt  =0,  -c  ^ x <_  c ; (3.24) 

' S2  a*: 


c 

/ A(s)KN3(c,y;s)ds 
-c 


c 

+ j B(s)K  (c,y;s)ds 

-c 


h _ 

+ / C(t) [LN3(c,y;t)  - LN3(c,y;t)]dt 


h _ 

+ / D(t)L  (c,y;t)dt  * 0,  0 < y < h i (3.25) 

0S3 


c 

/ AfsjK^C/yjsJds  + 
-c 


c 

/ B(s)K24(c,y;s)ds 
-c 


h _ 

+ / C(t) [-LN4(c,y;t) ]dt 
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+ / D(t)  [L  (c , y ; t)  + L (c,yjt) ]dt  = 0 0<_y<^h. 


(3.26) 


Note  that  Eqs . (3.6)- (3. 7)  are  eliminated  as  unnecessary.  It  is  also 
clear  that  A(x)  must  be  symmetric  and  B(x)  antisymmetric  in  the  range 
-c  <_  x <_  c.  These  properties  will  be  incorporated  in  the  numerical 
analysis . 

Since  the  unknown  stress  discontinuities,  A(s),  B(s),  C(t),  and 
D (t) , have  integrable  singularities  at  the  end  points,  in  addition  to 
Eqs.  (3. 23) - (3. 26) , they  must  also  satisfy  the  following  subsidiary 
conditions  [8] : 


c h 

/ A(x)dx  + 2 / D(y)dy  = 2F  ; 
-c  0 


(3.27) 


c 3u 

/ [-r*]  dx  = 0 

o 3*  y=° 


(3.28) 


c 3u  h 9u 

/ [-r-£]  . d x - / [-A  n dy  = 0 

h 3x  y=0  dy  x=0 


(3.29) 


B (0)  = 0 


(3.30) 


The  choice  of  Eqs.  ( 3. 27) - (3 . 30)  as  the  subsidiary  conditions 
appropriate  for  this  problem  warrants  a brief  discussion.  It  was 
found  that  extreme  care  must  be  exercised  in  selecting  the  correct 
combination  of  conditions  from  several  available  relationships  which 
could,  conceivably,  serve  as  subsidiary  conditions.  The  task  is  to 
prescribe  conditions  which  supplement  the  boundary  conditions  in 
describing  the  problem,  without  using  redundant  ones,  and,  of  course, 
without  violating  any  of  the  elasticity  rules  that  apply  to  mixed 
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boundary  value  problems.  Including  redundant  conditions  has  the  same 

effect  as  using  the  same  equation  twice  (i.e.  renders  the  matrix  of  ‘ 

coefficients  singular) , while  the  subtle  manner  in  which  certain  fac- 
tors enter  the  problem  creates  the  danger  of  violating  fundamental 
elasticity  rules,  such  as  overconstraining  the  problem  by  prescribing 
tractions  and  displacements  along  the  same  direction. 

The  governing  singular  integral  equations  (Eqs.  3.23  - 3.26) 
ensure  that  the  base  and  the  sides  of  the  block  do  not  rotate  and  do 

not  elongage  (shrink).  Equation  (3.27)  represents  the  global  equili-  ' 

brium  for  the  rigid  block  in  the  vertical  direction.  Equation  (3.30) 
is  a statement  requiring  B(x)  to  be  antisymmetric.  This,  along  with 
Eq.  (3.22)  which  requires  that  E(y)  = -C(y),  satisfies  the  equation 
of  equilibrium  in  the  horizontal  direction.  Equations  (3.28)  and 
(3.29)  impose  kinematic  restrictions  on  the  problem.  Equation  (3.28) 
fixes  the  sides  of  the  block  at  x = + c,  while  Eq.  (3.29)  requires  that 
the  base  and  the  sides  of  the  block  do  not  move  relative  to  each  other 
in  the  vertical  direction. 

Equations  (3. 23)- (3. 30)  are  therefore  the  relationships  appro- 
priate for  the  solution  of  the  problem  of  a rigid  block  embedded  in  an 
elastic  half  space  and  undergoing  vertical  displacement. 


Ai 


THE  TWO- INCLUSION  PROBLEM:  VERTICAL  TRANSLATION 


The  integral  equations  associated  with  the  vertical  line  inclu- 
sions embedded  in  an  elastic  half  space,  which  were  derived  in  Chapter 
II,  are  used  in  this  section  to  formulate  a problem  which  is  analogous 
to  the  rigid  block  one,  although  much  simpler.  This  is  the  problem  of 
two  inclusions  of  length  h,  parallel  to  each  other,  perpendicular  to 
the  free  surface  of  the  elastic  half  space  and  located  at  x *>  c and 
x = -c,  respectively.  Furthermore,  these  two  inclusions  are  thought 
of  as  being  rigidly  connected  to  each  other  above  the  surface,  so  that 
they  cannot  move  relative  to  each  other.  A load  P acting  in  the  nega- 
tive y-direction  is  applied  to  each  of  them,  and  they  are  to  translate 
vertically  without  rotation* 

The  boundary  conditions  appropriate  for  this  problem  are  given 


next : 


3u 

X 

3a11 

X 

3y 

3u 

. II 

3u 

Z 

9y  " 

3y 

3u 

X 

Bu11 

X 

3y  “ 

3y 

3u^ 

3uIX 

3u 


III 


(c,y)  + 


3y 

III 


(c,y ) = 0 


, 0 < y < h,  x = c ; (3. 31) 


9u 


(c,y)  + ~§y — (c,y)  = 0 


3u 


III 


(-c,y)  + 


9y 

III 


(-c,y)  * 0 


, 0<y<h,  x “ -c. 


3u 


(3.32) 


Equations  (3.31)  - (3.32)  represent  a system  of  four  singular 
integral  equations.  However,  the  symmetry  conditions  given  in  Eq. 
(3.22)  also  apply  to  the  two- inclusion  problem.  Therefore,  it  is  only 
necessary  to  solve  for  two  of  the  unknown  stress  discontinuities. 
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Equation  (3.22)  and  the  results  of  Chapter  II  are  incorporated  to 
write  Eqs.  (3.31)  in  the  following  form: 
h 

^ C(t> fLN3(C'y;t)  " LN3(c,y;t)Jdt 
h 

+ l°(t)  rs3(C'y,t)dt  “ 0.  0 < y < h, 


. •• 

/ C(t)(-L  (c,y;t)]dt  + / D(t)[L  (c,y,t)  + L (c,y:t) )dt  = 0, 

o o 

0 1 y 1 h-  (3.34) 

The  kernels  appearing  in  Eqs.  (3. 33)- (3. 34)  are  as  defined  earlier. 

Equations  (3. 33)- (3. 34)  may  be  solved  for  the  unknown  stress  dis- 
continuities, C (t)  and  D(t),  provided  that  the  subsidiary  conditions 
listed  next  are  satisfied: 

h 

/ D(y)dy  » P . (3. 35) 

0 ' 

c 3U 

/Hr-]  . dx  - 0.  (3.36) 

0 ax  y=0 

The  physical  significance  of  Eqs.  (3. 35)- (3. 36)  is  clear.  Equa- 
tion (3.35)  is  the  equation  of  equilibrium  in  the  y-direction,  while 
Eq.  (3.36)  fixes  the  inclusions  at  x • +c.  Note  that  the  equation  of 
equilibrium  in  the  x-direction  is  also  satisfied,  since  E(t)  =-C(t). 
Concerning  the  order  of  the  singularities  at  y * 0,h,  it  should  be 
noted  here  that  while  the  one  at  the  free  surface,  (+  c ,0),  remains 
unchanged  and  is  determined  from  Eq.  (3.14),  h = 0.71117,  the  singularity 
at  (+c,h)  appropriate  for  the  two-inclusion  problem  is  square  root. 


CHAPTER  IV 


STRESS  ANALYSIS  OF  EMBEDDED  BLOCKi  HORIZONTAL  DISPLACEMENT 

HORIZONTAL  TRANSLATION  OF  RIGID  BLOCK 

The  problem  investigated  in  this  section  is  that  of  a rigid 
block  partially  embedded  in  an  elastic  half  space  and  forced  to  trans- 
late in  the  negative  x-direction  without  rotation,  under  the  action  of 
a horizontal  load  2Q,  and  a moment  M,  as  shown  in  Fig.  3. 

The  boundary  conditions  appropriate  for  this  problem  are  the 
same  as  those  given  in  Chapter  III,  Eqs.  (3.1)— (3.3) , for  the  vertical 
translation  of  the  rigid  block  problem,  and  they  are  repeated  next: 

9u  9u 

°,  °.  -c  < x < c,  y = h;  (4.1) 

9u  9u 

-jj  = 0,  =0;  0 < y < h,  x = +c.  (4.2) 

Equations  (4.1)  - (4.2)  represent  a system  of  six  singular 
integral  equations  in  the  unknown  stress  discontinuities  A(s),  B (s ) , 

C (t) , ...,  F (t) , which  were  defined  in  Chapter  II.  The  symmetry  of 
the  problem  requires  that  the  following  relationships  hold: 

E(y)  = C (y) , F (y)  - -D(y),  (4.3) 

A (x)  = -A(-x),  B (x)  * B(-x).  (4.4) 

Equation  (4.3)  reduces  the  number  of  unknowns  from  six  to  four, 
thereby  allowing  the  elimination  of  two  of  the  governing  singular 
integral  equations.  Equation  (4.4)  will  be  incorporated  in  the 
numerical  analysis. 
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Using  Eq.  (4.3)  the  four  governing  equations  are  now  written 
in  the  following  form: 

c c 

/ A(s)K  (x,h;s)ds  + / B(s)K  (x,hjs)ds  + 

-c  Ni  -c  bL 


h 

+ / C(t) [L  . (x,h;t)  + I. (x,h;t)]dt  + 

q N 1 N 1 


h 

+ / D (t) [Lsl(x,h;t)  - Lgl(x,h;t) ]dt  = 0,  -c  £ x ^ c;  (4.5) 


c C 

/ A(s)KN2(x,h;s)ds  + / B(s)KS2 


-c 


-c 


(x,h;s)ds 


+ 


h _ 

+ / C(t) (LN2(x,h;t)  + LN2(x,h;t)]dt  + 


h _ 

+ / D(t)  [L  (x,h;t)  - Lc_  (x,  h;  t)  ]dt  «=  0,  -c  £ x < c;  (4.6) 

0 2 2 


c c 

/ A(s)K_(c,y;s)ds  + / B (s )Kg3 (c ,y ; s ) ds  + 
-c  -c 

h _ 

+ / C(t) [LN3(c,y;t)  + LN3(c,y>t) ]dt  + 


h _ 

+ / D (t)  [~L„  (c,y ; t ) ] dt  = 0,  0 <_  y £ h>  (4.7) 


r 


: 

! .. 

* 

! 

I 
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c c 

/ A(s)KN4 (c,y;s)ds  + / B(s)Kg4 


-c 


-c 


(c,y > s) ds 


+ 


h _ 

+ / c(t)L  (c,y;t)dt  + 

0 

h _ 

+ / D (t) [L  (c,y ; t)  - L 4(c,y;t)]dt  = 0,  0 £ y £ h.  (4.8) 

0 S'*  S'* 

The  kernels  K K_  . , L , and  L.,,  i = 1,  . ..,  4,  are  given  in 
Ni  oi  Ni  si 

Appendix  A.  The  barred  kernels  L .,  L , i=l,  4,  are  obtained 

Ni  Si 

from  the  corresponding  unbarred  kernels  by  replacing  c by  -c. 

Note  that  the  order  of  the  singularities  at  the  surface,  (+c,0), 
and  at  the  corners,  (+c,h) , remains  unchanged  and  is  as  given  in 
Chapter  III,  Eqs.  (3.14)  and  (3.21).  Of  course,  the  same  holds  true 
for  the  corner  conditions,  and  the  ones  appropriate  for  the  corner 
(c,h)  are  given  by  Eq.  (3.16). 

Equations  (4.5)- (4.8)  can  be  solved  for  the  unknown  stress  dis- 
continuities A(s),  B(s),  C (t) , and  D(t),  provided  that  a set  of  sub- 
sidiary conditions  is  also  satisfied.  The  need  for  additional  condi- 
tions is  due  to  the  fact  that  A(s),  ...,  D(t)  have  integrable  singu- 
larities at  the  end  points,  as  explained  in  Chapter  III.  The  sub- 
sidiary conditions  associated  with  Eqs.  (4. 5) -(4.8)  are  given  nexti 

c h 

/ B (x ) dx  + 2/c(y)dy  - 2Q  ; (4.9) 

-c  0 


c 3u 

/ >-# 


x'y-0 


dx  = 0 ; 


(4.10) 
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C du 

I l X 


0 


A (0) 


x y=0 


0 . 


h 8u 

dx-  / [ X 

0 


“57]x=o  dy  = 0 ' 


(4.11) 


(4.12) 


i 


The  discussion  presented  in  Chapter  III  concerning  the  subsi- 
diary conditions  (3. 27)- (3. 30) , also  applies  to  the  selection  of  Eqs. 
(4.9)-(4.12)  as  the  correct  set  of  conditions  needed  to  supplement  the 
governing  singular  integral  equations,  Eqs.  (4.5)- (4.8),  and  will  not 
be  repeated  here. 

It  suffices  to  state  that  Eqs.  (4.9)  and  (4.12),  the  latter  in 
association  with  Eq.  (4.3),  represent  the  global  equations  of  equili- 
brium for  the  block  in  the  x-  and  the  y-  directions,  respectively. 
Equation  (4.10)  ensures  that  the  sides  of  the  block  do  not  translate 
vertically,  while  Eq.  (4.11)  requires  that  the  base  and  the  sides  of 
the  block  do  not  move  relative  to  each  other  in  the  horizontal  direction. 


THE  TWO-INCLUSION  PROBLEM : HORIZONTAL  TRANSLATION 


In  Chapter  III  the  two-inclusion  problem  was  described,  and  for- 
mulated for  the  vertical  translation  case.  The  present  section  con- 
siders the  two- inclusion  problem  undergoing  horizontal  translation, 
under  the  action  of  a load  2Q  in  the  negative  x-direction.  As  pre- 
viously, the  rigid  inclusions  are  thought  of  as  being  rigidly  connec- 
ted to  each  other  above  the  surface,  and  they  are  to  translate  without 
rotation. 

The  boundary  conditions  for  this  problem  are  identical  to  the 
ones  given  by  Eqs.  (3. 31)- (3. 32) , and  are  repeated  next: 
du  3u 

~ 5y  “ ~~$y  * 0 £ y 5.  x = ±c  • (4.13) 

The  symmetry  relationships  given  in  Eq.  (4.3)  are  incorporated 
to  write  the  governing  singular  integral  equations,  at  x = c,  in  the 
following  form: 
h 

^ C(t) (LN3(c,y;t)  + LN3(c,y;t) ]dt  + 
h _ 

+ / D(t) [-L_  (c,y;t)]dt  * 0,  0 < y < h;  (4.14) 

0 S3  " ~ 

h 

/ C(t)  L (c,y;t)dt  + 

0 

h _ 

+ / D(t)  [L  (c,y;t)  - L . (c,y?t)  ]dt  * 0,  0 <_  y <_  h.  (4.15) 

0 5,4 


23 


24 


The  stress  discontinuities  C(t)  and  D(t)  were  defined  in  Chapter  II, 

Eq.  (2.11),  while  the  kernels  ^ and  Lg^  ^ are  given  in  Appendix  A. 
Recall  that  the  barred  kernels  are  obtained  from  the  corresponding 
unbarred  ones  by  replacing  c by  -c. 

The  functions  being  sought,  C(t)  and  D(t),  have  integrable  sin- 
gularities at  the  end  points;  therefore,  besides  equations  (4.14)- (4.15) 
they  must  also  satisfy  the  following  subsidiary  conditions: 

h 

£ C(y)dy  = Q ; (4.16) 

c 3u 

/ l-£l  0 - °-  (4-17> 

0 

Equation  (4.16)  is  the  equation  of  equilibrium  in  the  x-direc- 
tion,  while  equilibrium  in  the  y-direction  is  also  satisfied  by  vir- 
tue of  the  symmetry  of  the  problem,  Eq.  (4.3).  Equation  (4.17)  requires 
that  the  rigid  inclusions  do  not  translate  in  the  y-direction. 

In  regard  to  the  order  of  the  singularities  for  this  problem, 
note  that  the  singularity  at  (+c,h)  is  square  root,  while  that  at  the 
surface,  (+c,0),  remains  unchanged,  n = 0.71117. 


CHAPTER  V 


STRESS  ANALYSIS  OF  EMBEDDED  BLOCK:  ROTATION 


ROTATION  OF  RIGID  BLOCK 


This  chapter  investigates  the  problem  of  a partially  embedded 
rigid  block  that  rotates  through  an  angle  a under  the  action  of  a 
moment  M,  as  shown  in  Fig.  3.  Here,  the  tangential  load  Q is  set 
equal  to  zero.  The  boundary  conditions  for  this  problem  are  given  next: 

3u  3u 

~Am  °'  °»  -c  -x-c'  y = h'  (5-1} 

3u  3u 


— ~ ” -a#  • 0;  0 < y < h,  x = +c.  (5.2) 

3y  3y  — — — 

Substitution  of  the  horizontal  and  vertical  inclusion  results 
derived  in  Chapter  II  in  Eqs.  (5.1)- (5.2)  leads  to  a system  of  six 
governing  singular  integral  equations  with  generalized  Cauchy  kernels 
However,  the  symmetric  nature  of  the  problem  allows  for  the  elimina- 
tion of  two  of  the  unknown  stress  discontinuities,  since  the  relation 
ships  given  next  must  be  satisfied: 

E(y)  - C(y);  F(y)  - -D(y).  (5.3) 

Therefore,  only  four  of  the  equations  represented  by  Eqs.  (5.1) -(5. 2) 
are  necessary  for  the  solution  to  this  problem.  Note  that  the  sym- 
metry of  the  problem  also  requires  that  A(x)  is  antisymmetric,  and 
B (x)  symmetric  in  the  range  -c  < x < c. 
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Equation  (5.3)  is  incorporated  to  write  the  equations  represented 
by  Eq.  (5.1)  and  those  of  (5.2)  at  x = c in  the  following  form: 


4x(~+l)7  { 7 MsjK^tx.hjsJds  + / B(s)l^1(x,h;S)ds 


h _ 

+ / C(t) [LN1(x,hjt)  + LN1(x'h;t) Idt 


h 

+ / D (t) [L  (x,h jt)  - L (x,hj t) ]dt}  = 0,  -c  < x < c:  (5.4) 

0 S1  Si 


- { / A (s) K (x,h;s)ds  + / B(s)K  (x,h;s)ds 

4tt(k+1)v  N2  ' S l 


h _ 

+ / C(t) (LN2(x,h;t)  + LN2(x,h;t) ]dt 


h _ 

+ / D (t)  [L  _ (x , h; t)  - L (x,h>t)]dt)  = a,  -c  £ x _<  c;  (5.5) 


C C 

4^(7^.!)-  <J  A(s)KN3(c,y,s)ds  + Jb  (s)  Kg3  (c,y;  s)  ds 


h _ 

+ / C(t)  [LN3(c,y;t)  + L^^yjtHdt 


h 

+ / D(t) l-L  (c,y;t)]dt}  * -a,  0 <y  < hi 

0 SJ 


(5.6) 
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4ir  (ic+1) p J A(s)KN4(c,y;s)ds  + / B (s)  (c,y , s) 
“ c * c 


h 

+ / c(t)L  (c,y;t)dt  + 
0 N 


h 

+ / D(t) [Lg4(c,y;t)  - Lg4 (c,y; t) ] dt}  =0,  0 < y < h.  (5.7) 

The  kernels  appearing  in  Eqe.  (5.4)— (5. 7)  are  as  defined  previously. 
Note  that  the  corner  conditions  for  this  problem  are  given  by  Eqs. 
(3.16),  for  the  corner  (c,h) . The  order  of  the  singularities  at  the 
points  (+c,0),  and  (+c,h)  was  derived  previously,  and  is  given  by 
Eqs.  (3.14)  and  (3.21),  respectively. 

The  unknown  functions  A(s),B(s),  C(t),  and  D(t)  have  integrable 
singularities  at  the  end  points.  Thus,  Eqs.  (5. 4) -(5. 7)  must  be  sup- 
plemented by  additional  conditions,  and  the  ones  appropriate  for  this 
problem  are  given  next: 

c h 

/ B (x)dx  + 2 / C(y)  dy  « 0 I (5.8) 

-c  0 


c 3u 

/ dx  = ca  ; (5.9) 

0 3x  y=0 


l h£’y-0  0*  - / d>  ■ h“  1 <5'10’ 

A (0)  - 0.  (5.11) 

Equation  (5.8)  represents  the  global  equilibrium  equation  for  the 


rigid  block  in  the  x-direction.  Equilibrium  in  the  y-direction  is 
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satisfied  by  Eq.  (5.11),  which  requires  that  A(x)  is  antisymmetric, 
in  association  with  Eq.  (5.3),  requiring  that  F(y)  = -D(y).  Equa- 
tions (5.9)-(5.10)  impose  kinematic  restrictions  on  the  problem. 
Equation  (5.9)  prescribes  the  translation  of  the  corner  (c,0)  of 
the  block  in  the  y-direction,  while  Eq.  (5.10)  gives  the  relationship 
between  the  horizontal  displacements  of  points  (c,0)  and  (0,h)  on  the 
block.  It  should  be  noted  that  the  usual  small  angle  approximations 
have  been  employed  in  writing  Eqs.  (5.9) - (5. 10) . 


THE  TWO- INCLUSION  PROBLEM;  ROTATION 

The  problem  considered  in  this  section  is  that  of  two  vertical 
rigid  inclusions,  partially  embedded  in  an  elastic  half  space  and  per- 
fectly bonded  to  it,  as  described  earlier  in  Chapter  III,  and  loaded 
in  such  a manner,  so  that  they  both  rotate  through  an  angle  a. 

The  boundary  conditions  for  this  problem  are  given  next: 

9u  9u 

— — = -a,  — =0;  0 < y < h,  x = +c.  (5.12) 

9y  3y  _ _ — 

The  symmetry  conditions  given  in  Eq.  (5.3)  reduce  the  number  of  un- 
known stress  discontinuities  to  two;  hence,  only  two  of  the  four 
singular  integral  equations  represented  by  Eq.  (5.12)  are  needed  for 
a solution. 

Equation  (5.3)  is  used  to  write  the  equations  governing  this 
problem,  at  x = c,  in  the  following  form: 
h 

h 

+ ^ D(t)  l-Ls3(c,y;t)  ]dt }=*-<*,  0 f.  y £ h;  (5.13) 

h 

< £ =<«  vc'”t,at 

h 

+ ^ D (t)  [LS4(c,y;t)  - Ls4(c,y;t)]dt}  = 0,  0 <_  y <_  h.  (5.14) 

Equations  (5. 13)- (5. 14)  may  be  solved  for  the  unknown  functions 
C(t)  and  D(t),  provided  that  the  subsidiary  conditions  given  next  are 


satisfied: 
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h 

/ C(y)dy  = 0 ; (5.15) 

0 

c 3u 

/ [— r^-]  dx  = ca  . (5.16) 

0 3x  y=  0 

Equation  (5.15)  is  the  equation  of  equilibrium  in  the  x-direction, 
while  Eq.  (5.16)  is  identical  to  Eq.  (5.9),  which  was  written  for  the 
rigid  block  problem,  and  it  prescribes  the  displacement  in  the  y- 
direction  of  the  point  (c,0)  of  the  block  relative  to  the  midpoint 
at  the  surface. 

As  in  the  previous  cases  of  the  two-inclusion  problem,  the  singu- 
larities are  of  order  (n-1)  at  the  surface  points  (+c,0),  and  (0-5-1) 
at  the  points  (+c,h) , where  n = 0.71117. 


A4 


CHAPTER  VI 


NUMERICAL  ANALYSIS  AND  DISCUSSION 

NUMERICAL  ANALYSIS 

The  objective  of  the  numerical  analysis  is  to  solve  numerically 
the  systems  of  governing  singular  integral  equations  with  the  corres- 
ponding subsidiary  conditions,  derived  in  the  previous  three  chapters, 
for  the  unknown  stress  discontinuities. 

The  equations  are  normalized  by  introducing  the  following  variable 
changes : 

x = cx,  s = cs;  (6.1) 

y = |d+y)  , t = |(l+t)  . (6.2) 

In  addition,  the  stress  discontinuities  are  given  forms  that  reflect 
their  correct  singularities  at  the  corners  (and  the  surface) . For  the 
case  of  vertical  displacement  of  the  rigid  block  the  following  substi- 
tutions are  made: 

A (s ) ■=  - A(s)  (l-I2)  , B (s ) = ^B(s)  (1-72)  ; (6.3) 

c c 

c (t)  = ^ C(t)  d-t)  (l+t)0"1? 

2p -c-i  — n_!  <6'4> 

D (t)  = D(t) (1-t)  ^A(l+t)n  . 
n 

Equations  (6.1)-(6.4)  also  apply  to  the  horizontal  translation  of 
the  rigid  block  problem,  provided  that  P is  replaced  by  Q.  Note  that 
Eqs.  (6.2)  and  (6.4)  are  also  employed  to  normalize  the  equations 
derived  for  the  two-inclusion  problem,  for  the  vertical  and  horizontal 


translation  cases 
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To  normalize  the  equations  derived  for  the  case  of  rotation  of 
the  rigid  block  (see  Chapter  V)  the  following  substitutions  are  needed: 


A (s) 

B (s) 

C(t) 


2u  (K+l)yac  r.~  -2.C-1 

— - A (s)  (1-s  ) ; 


2n  (x+1) pac  -2. C-l 

B (s ) (1-s  ) ; 


c(t)  (i-t)c"1(i+t)n’1l 


(6.5) 


(6.6) 

o(t)  - 5(I)(i-t){-1(itt)H  . 

K 

The  substitutions  given  in  Eqs.  (6.1)-(6.2)  also  apply  to  this  problem. 
Equations  (6.2)  and  (6.6)  are  also  used  to  normalize  the  equations 
associated  with  the  two- inclusion  problem  undergoing  rotation. 

Equations  (6.1) -(6. 6)  will  appropriately  normalize  the  governing 
singular  integral  equations  and  the  subsidiary  conditions  so  that  the 
numerical  analysis  may  be  conducted.  An  adjustment  is  required  in 
Eqs.  (6.1)  for  normalizing  Eqs.  ( 3. 28) - (3. 29) , (3.36),  (4. 10) - (4. 11) , 
(4.17),  (5.9)- (5.10) , and  (5.16);  instead  of  x = cx,  take 


x = c(l+x)/2.  (6.7) 

For  convenience,  the  parameter  y is  introduced,  and  is  defined  as 
Y “ c/h.  (6.8) 

It  should  be  noted  that  for  the  case  of  rotation  of  the  rigid  block, 
it  becomes  necessary  to  assign  a specific  numerical  value  to  either 
c or  h.  Thus,  for  that  problem  h is  given  unit  length,  and  therefore 
Y becomes  equal  to  c. 

The  numerical  scheme  to  be  used  in  this  analysis  is  that  described 
in  Erdogan,  Gupta,  and  Cook  [8] , and  is  a collocation  scheme  based 
upon  formulas  for  Gaussian  Integration  (see  Stroud  and  Secrest  [10]). 
The  integration  points,  s^,  t^»  and  the  collocation  points,  x^,  y^, 
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are  defined  as  the  roots  of  the  Jacobi  polynomials  as  follows: 

(6.9) 


- 0;  - 0, 


PB«— 


(6.10) 


The  corresponding  weights  are  obtained  from  the  Gauss-Jacobi  integra- 
tion formula  (see  [10]). 

Note  that  since  the  subsidiary  conditions  which  impose  kinematic 
restrictions  on  all  problems  require  integrations  to  be  performed 
along  y = 0 and  x = 0,  additional  integration  points  and  weights  are 
needed.  The  integration  points  appropriate  to  the  interval  (0,c) 


along  y=*0  are  obtained  from 

„ (n-1,0)  - 


N 


(X.)  = o. 


(6.11) 


while  the  ones  for  (0,h)  along  x * 0 are  obtained  from  the  Legendre 
polynomials,  as  the  Gauss  formula  for  numerical  integration  may  be 
used  in  this  case,  since  there  are  no  singularities  at  the  end  points 
of  the  interval.  The  zeros  of  the  Legendre  polynomials  (integration 
points),  and  the  weights  associated  with  the  Gauss  formula  can  be 
obtained  from  several  sources  (see  e.g.  Abramowitz  and  Stegun  [11]). 

The  corner  conditions  are  in  the  forms  given  by  Eqs.  (3.16)- 
(3.17),  where,  for  the  cases  of  vertical  and  horizontal  translation, 
the  following  relationships  are  valid: 


Aq  - A(1)/2V,  Bq  = B(1)/2V; 

Aq  - A(-1)/2V,  - B(-1)/2V; 

WVFo " ^(1)'  °(1)'  ®(1>'  *’(1)- 


(6.12) 

(6.13) 


For  the  rotation  of  the  rigid  block  problem  Eq.  (6.12)  is  replaced  by: 
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RESULTS  AND  DISCUSSION 


The  numerical  analysis  was  completed  for  a range  of  the  geo- 
metric parameter  y (0.05  <_  y <_  8)  , and  for  Poisson's  ratio  (plane 
strain)  of  v = 0.3.  The  rate  of  convergence  was  tested  by  varying 
the  number  of  points  used  at  eight  point  intervals,  from  N = 16  to 
N = 48.  The  convergence  of  the  global  results  appeared  satisfactory 
and  will  be  discussed  in  detail  later.  All  results  presented  in 
this  dissertation  will  be  given  for  N = 48. 

The  results  for  the  respective  problems  of  vertical  displacement, 
horizontal  translation  without  rotation,  and  rotation  of  an  embed- 
ded rectangular  block  will  be  given  next. 

Vertical  displacement.  The  load  diffusion  curves  for  several 
values  of  y are  shown  in  Fig.  4 for  the  case  when  the  block  moves  ver- 
tically without  rotation.  The  load  acting  on  the  block  at  a distance 
y below  the  surface  is  given  as  a fraction  of  the  total  load,  where 
both  load  and  distance  have  been  put  in  dimensionless  form.  The  load 
given  by  the  intersection  of  the  load-diffusion  curves  and  the  y = h 
axis  represents  the  load  carried  by  the  base  of  the  block  for  each  y. 
Analogous  to  the  block  problem  is  that  of  two  parallel  inclusions, 
each  subjected  to  a vertical  load  P and  having  no  rotation.  The 
results  for  this  problem  can  be  compared  with  those  of  the  vertical 
displacement  for  the  block.  The  load  diffusion,  calculated  for  the 
two- inclusion  case,  is  shown  in  Fig.  5,  where  it  is  compared  to  the 
rigid  block  results;  the  solid  lines  represent  values  obtained  from 
Fig.  4,  and  the  dashed  lines  represent  the  two-inclusion  case.  These 
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latter  curves  go  to  zero  at  y = h,  since  there  is  no  base  to  carry  part 
of  the  load.  Furthermore,  the  dashed  and  solid  curves  are  relatively 
close  provided  that  0 < y < 0.8h.  Thus,  for  relatively  slender  blocks, 
a solution  which  ignored  the  base  would  give  approximately  valid  re- 
sults. Figure  6 shows  the  percentage  of  the  applied  load  P carried 
by  the  side  of  the  rectangular  block  as  a function  of  c/h.  As  ex- 
pected, when  c/h  is  small  most  of  the  load  is  carried  by  the  sides;  as 
c/h  increases  the  base  carries  a larger  share  of  the  load.  When  Y is 
approximately  equal  to  0.42,  the  sides  and  the  base  each  carry  half 
the  load. 

Figures  7 and  8 show  the  variations,  respectively,  of  the  discon- 
tinuities in  the  normal  stresses  along  the  base  of  the  block,  and  those 
in  the  shear  stresses  along  the  sides  for  several  values  of  c/h.  Since 
a zero  stress  state  prevails  throughout  the  rigid  block,  these  stress 
discontinuities  represent  the  actual  stress  distributions  along  the 
walls  of  the  elastic  body  surrounding  the  block,  as  is  evident  from 
the  definitions  of  A(x),  B(x),  C(y),  and  D(y).  Also,  since  the  normal 
stresses  are  symmetric  along  the  base,  they  are  plotted  only  for 
0 < x < c and  exhibit  singular  behavior  at  the  end  x = c.  The  shear 
stresses  are  plotted  for  0 < y < h,  and  singular  behavior  is  noted 
near  both  ends,  y = 0,h.  It  is  further  observed  that  the  magnitudes 
of  both  A (x)  and  D(y)  vary  proportionately  with  the  amount  of  load 
carried  by  the  base  and  the  sides,  respectively.  Thus,  the  magnitude 
of  A (x)  increases  with  y,  while  the  magnitude  of  D(y)  decreases  as  y 
increases. 


One  estimate  of  the  accuracy  of  the  results  may  be  obtained  by 
calculating  the  displacements  in  the  y-direction  along  the  surface  of 
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the  block  y = 0,  relative  to  the  midpoint  (0,0).  Since  the  block 
translates  as  a rigid  body  along  y,  these  displacements  should  be 
small  when  compared  to  those  obtained  from  the  solution  of  the  two- 
inclusion  problem  undergoing  vertical  displacement.  The  results  of 
this  comparison  are  given  in  Table  1.  It  is  noted  that  the  displace- 
ments at  the  corner  x = c for  the  two  inclusion  problem  vary  from 
being  approximately  18  times  greater  than  the  block  displacements 
(for  y * 0.3),  to  being  175  times  greater  (for  y = 4) . The  dis- 
placements obtained  for  the  rigid  block  problem  tended  to  oscillate 
about  zero,  while  the  ones  for  the  two  inclusion  problem  grew  with 
distance  from  (0,0),  as  expected.  It  should  be  added  that  for 

y < 0.3,  the  value  of  the  ratio  |u  /u  | approaches  unity.  That  is, 

yb  y2 

as  the  block  becomes  slender,  the  presence  (or  absence)  of  the  base 
has  little  or  no  effect  on  the  results  calculated  near  the  surface. 
Therefore,  a comparison  of  those  results  for  the  purposes  of  Table  1 
would  have  little  meaning,  since  the  surface  displacements  are  already 
small  to  within  the  accuracy  of  the  solution. 

Horizontal  displacement  without  rotation.  Figures  9-12  present 
results  obtained  from  the  solution  of  the  embedded  block  undergoing 
horizontal  displacement.  Figure  9 shows  the  variation  with  y of  the 
percentage  of  the  total  load  that  is  carried  by  the  side  of  the  rigid 
block.  In  this  case,  the  point  at  which  half  of  the  load  is  carried 
by  the  sides  and  half  by  the  base  is  approximately  Y = 1.5. 

The  shear  stress  distribution  along  the  base  of  the  block,  and 
the  normal  stress  distribution  along  the  side  x = c of  the  block  are 
shown  for  several  values  of  Y in  Figs.  10  and  11,  respectively.  The 
results  showing  the  variation  of  the  magnitudes  of  the  stresses  with 
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Y are  analogous  to  the  results  for  the  vertical  displacement  case:  the 
magnitude  of  the  stresses  tends  to  increase  with  the  amount  of  load 
carried  by  the  surface  of  the  block  upon  which  they  act. 

Since  the  block  translates  horizontally  without  rotation,  an  ap- 
plied moment,  M,  in  addition  to  a horizontal  load,  2Q,  must  be  ap- 
plied to  prevent  the  block  from  rotating.  The  ratio  M/Qh  is  plotted 
in  Fig.  12  as  a function  of  c/h.  At  the  limit,  as  c/h  approaches  zero, 
the  results  are  in  agreement  with  the  result  obtained  from  the  solu- 
tion for  a single  inclusion  perpendicular  to  the  free  surface,  y = 0, 
and  loaded  by  a horizontal  force,  2Q,  and  a moment,  M,  such  that  the 
inclusion  translates  in  the  negative  x-direction  without  rotation. 

Rotation.  Results  obtained  for  the  problem  of  rotation  of  the 
embedded  rigid  block  are  presented  in  Figs.  13-17.  The  normal  and 
shear  stresses  acting  on  the  base  of  the  block  are  given  in  Figs.  13- 
14  as  functions  of  the  distance  from  the  midpoint,  for  several  values 
of  c/h.  The  normal  and  shear  stress  distributions  along  the  sides  of 
the  block  for  0 < y < h,  are  plotted,  respectively,  in  Figs.  15-16.  The 
curves  show  the  variation  of  the  magnitude  and  sign  of  the  stresses 
with  y.  A better  understanding  of  the  mechanics  of  the  problem  can  be 
obtained  by  examining  the  moment  equilibrium  equation  for  the  block, 
which  is  given  next: 

c c h 

M = -/  A (x)x  dx  + h / B (x)  dx  + 2 / C(y)y  dy 

-c  -c  0 

h 

-2c  | D (y)  dy.  (6.15) 

0 

In  writing  Eq.  (6.15)  use  has  been  made  of  Eq.  (5.3).  Using  the  sym- 
metries of  A (x)  and  B(x)  appropriate  for  the  rotation  problem,  and 
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incorporating  Eqs.  (6.1)-(6.2),  and  (6.5)-(6.6),  Eq.  (6.15)  is 

written  in  the  normalized  form  given  next: 

1 M tt(ic+1)  „ . 3 r T,-  ~ ,,  -2.5-1.- 
— j ~ = — [-4y  J A(x)  x (1-x  ) dx 


lih 


+ 4y2  / B (x) (1-x2) C-1dx 
0 

+ / C(7)  (1+y)  (l-y)C_1(H-y)n~1d7 
-1 

1 


-2t  / D(y)(l-y)5~  (1+y) ' 

-1 


dy] 


(6.16) 


Equation  (6.16)  may  be  written  as 


— — = a + B + C + D. 

, 2 a S S S S' 
yh 


(6.17) 


where  the  definitions  of  the  quantities  A , . ..,  D are  readily  ob- 

b b 

2 

tained  by  comparing  Eqs.  (6.16)  and  (6.17).  The  ratio  M/ph  a repre- 
sents a measure  of  the  rotational  stiffness  of  the  partially  embedded 
rigid  block.  Therefore,  the  quantities  A , Bc,  C , and  P represent 
the  contribution  made  to  the  stiffness  by  the  stresses  A(x),  B(x), 

C(y),  and  D(y),  respectively.  Table  2 shows  the  variation  of  these 
quantities  with  y.  Note  that  the  contribution  of  the  normal  stresses 
acting  on  the  base  of  the  block  becomes  significant  only  for  y > 1, 
while  the  opposite  is  true  for  the  normal  stresses  acting  on  the 
sides  (their  contribution  becomes  more  important  as  y becomes  smaller). 
It  is  also  observed  that  the  contribution  of  the  base  shear  stresses 
to  the  total  stiffness  remains  proportionately  low  for  all  values  of 

y,  while  that  of  the  side  shear  stresses  increases  with  y. 

2 2 2 
By  dividing  M/ph  a by  y , one  obtains  a new  ratio,  namely  M/pc  a. 

These  two  ratios  are  plotted  as  functions  of  y in  Fig.  17.  There  are 
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two  limiting  tests  for  these  results:  1)  as  c/h  tends  to  zero,  the 
2 

results  for  M/ph  a tend  to  the  single  inclusion  result;  2)  as  c/h 
becomes  large  the  results  may  be  compared  to  the  one  obtained  by 
Muskhelishvili  for  the  rotation  of  a rigid  stamp  bonded  to  an  elastic 
body  (see  [13],  p.  492,  Eq.  114.19a).  As  may  be  seen  in  Fig.  17, 
the  two  limiting  cases  seem  to  support  the  accuracy  of  the  calculated 
results. 

The  corner  conditions  derived  for  the  corner  (c,h)  were  checked 
by  calculating  the  values  of  the  stress  discontinuities  at  that 
corner,  and  substituting  in  Eqs.  (3.16).  The  corner  values  were  cal- 
culated by  a quadratic  extrapolation  applied  to  the  three  points 
nearest  that  corner.  It  was  found  that  the  corner  conditions  were  not 
satisfied;  the  size  of  the  error  was  generally  of  the  same  order  of 
magnitude  as  some  of  the  individual  terms.  Even  though  an  allowance 
must  be  made  for  the  error  introduced  by  the  quadratic  extrapolation, 
this  size  error  is  still  considered  large  for  the  present  type  of 
analysis.  To  correct  this  discrepancy,  and  more  importantly,  to  check 
the  accuracy  of  the  global  results,  three  of  the  corner  conditions 
given  in  Eqs.  (3.16)  were  incorporated  in  the  system  of  3N  equations, 
by  removing  the  equations  closest  to  the  corner  (c,h)  from  the  cor- 
responding boundary  conditions.  This  new  system  of  equations  yielded 
results  which,  away  from  the  corner  (c,h)  matched  the  results  pre- 
viously obtained  very  closely  (to  four  significant  figures) , and  in 
addition,  satisfied  all  the  corner  conditions.  Moreover,  the  global 
results  obtained  using  this  new  system  of  equations  varied  slightly 
from  the  ones  obtained  previously,  as  it  will  be  shown  later.  To 
further  test  the  sensitivity  of  the  global  results  to  the  collocation 


41 


scheme  used,  the  order  of  the  singularity  at  (c,h)  was  changed  from 
C - 1 = -0.40484  to  -0.5,  and  the  calculations  were  repeated  with 
N = 48.  Lastly,  the  singularities  at  (c,h)  and  (c,0)  were  set  equal 
to  -0.5,  and  a new  set  of  results  obtained  for  N = 48.  A comparison 
of  some  of  the  results  obtained  from  these  solutions  is  given  in 
Tables  3-4.  An  examination  of  these  tables  shows  that  the  global 
results  remain  relatively  unaffected  by  small  changes  in  the  order  of 
the  corner  singularities. 

It  should  be  noted  that  when  the  order  of  the  singularities  at 
both  (c,h)  and  (c, 0)  was  taken  as  -0.5,  the  related  Jacobi  polynomials 
reduced  to  the  Chebyshev  polynomials  of  the  first  kind  (see  e.g.  [11]). 
Numerical  solutions  of  Cauchy-type  singular  integral  equations  with 
regular  kernels  obtained  by  use  of  the  Gauss-Chebyshev  integration 
formula,  have  been  shown  to  converge  to  the  correct  results  by  Erdogan 
and  Gupta  [14] , and  Kalandiya  [15] . Furthermore,  Kalandiya  states  in 
[15]  that  this  method  was  also  successfully  applied  to  problems  with 
generalized  kernels.  Of  course,  as  the  singularities  used  here  are 
not  of  the  correct  order,  this  method  cannot  be  expected  to  yield  cor- 
rect results  at  the  corners.  However,  the  fact  that  the  global  re- 
sults obtained  by  use  of  the  Gauss-Chebyshev  integration  formula  match 
closely  those  obtained  using  the  Gauss-Jacobi  integration  formula,  as 
may  be  seen  in  Tables  3-4,  is  supportive  of  the  validity  of  the 
global  results. 

Although  the  accuracy  of  the  global  results  obtained  by  use  of 
the  collocation  scheme  introduced  in  [8]  has  been  shown  to  be  satis- 
factory, the  accuracy  of  the  results  in  the  vicinity  of  the  corners  is 
not  certain.  There  are  two  reasons  for  this:  1)  A second  root  is 
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present  at  the  corner  (c,h)  and  is  not  taken  into  account;  this  root 
has  a singular  contribution  of  order  0.75904-1  = -0.24096  near  that 
corner.  As  pointed  out  by  Westmann  in  [7],  although  the  stress  field 
in  the  vicinity  of  the  corner  is  dominated  by  the  largest  singularity, 
the  presence  of  two  singular  terms  may  be  important  to  the  problem. 

2)  In  two  recent  papers,  [16,  17],  Theocaris  and  Ioakimidis  suggest 
that  the  collocation  points  as  determined  for  the  Gauss-Jacobi  method 
may  not  be  the  correct  ones.  In  [16],  it  is  suggested  that  the  rate 
of  convergence  of  the  results  can  be  improved  by  a different  choice  of 
collocation  points  associated  with  the  Gauss-Jacobi  integration  rule, 
or  by  using  a different  scheme  such  as  the  Lobatto  rule.  The  colloca- 
tion points  for  the  case  of  real  singularities  are,  respectively,  the 
roots  of  the  Jacobi-functions,  and  of  their  derivatives  (see  Elliot 
[18]).  In  [17],  the  problem  of  an  antiplane  shear  crack  terminating 
at  a bimaterial  interface  solved  by  Erdogan  and  Cook  [19] , was  recon- 
sidered, and  numerical  results  were  obtained  by  using  the  modified 
Gauss-Jacobi  and  the  Lobatto-Jacobi  methods  discussed  in  [16].  It 
was  found  that  while  in  general  the  results  of  [17]  matched  closely 
those  of  [19],  the  value  of  the  stress  intensity  factor  at  the  inter- 
face obtained  by  Theocaris  and  Ioakimidis  was  in  much  better  agree- 
ment with  the  theoretical  result  than  that  obtained  by  Erdogan  and  Cook. 
However,  the  rate  of  convergence  was  not  satisfactory.  The  authors 
attributed  this  to  the  existence  of  a second  singularity  of  positive 
order  near  the  interface.  A technique  is  proposed  for  accounting 
for  both  singularities  in  the  analysis,  and  this  results  in  a much 
faster  convergence  of  the  results  to  the  theoretical  value.  However, 
it  should  be  realized  that  this  technique  could  not  be  applied  to  the 
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present  analysis,  since  it  deals  with  singularities  of  opposite  signs. 

In  view  of  the  arguments  just  presented,  the  accuracy  of  the  cor- 
ner values  of  the  stress  discontinuities  obtained  in  the  present  ana- 
lysis, by  use  of  the  collocation  scheme  introduced  in  [8],  is  uncer- 
tain. However,  for  completeness,  a discussion  of  these  values  is  pre- 
sented in  Appendix  B. 

To  investigate  the  rate  of  convergence  of  the  global  results,  the 

number  of  points  taken  was  varied  at  eight  point  intervals,  from 

N=16  to  N=48 . Several  of  the  calculated  global  results  are  presented 

in  Tables  5-7.  Tables  5-6  show  the  calculated  values  for  the  portion 

of  the  load  carried  by  the  sides  of  the  block  for  the  vertical  and 

horizontal  displacement  cases,  respectively.  The  variation  of  the 
2 

ratio  M/ph  a with  N,  is  shown  for  several  values  of  c/h  in  Table  7. 

It  is  seen  that  the  rate  of  convergence  of  these  representative  global 
quantities  is  quite  satisfactory. 


CHAPTER  VII 


CONCLUSIONS 

The  formulation  of  the  boundary  value  problem  of  a rigid  rec- 
tangular block  partially  embedded  in,  and  perfectly  bonded  to  an 
elastic  half  space  led  to  a system  of  singular  integral  equations, 
whose  unknowns  were  the  discontinuities  in  the  stresses  between  the 

bonded  surfaces.  Three  separate  loading  conditions  were  applied  , 

causing  the  block  to  1)  translate  vertically  without  rotation, 

2)  translate  horizontally  without  rotation,  and  3)  rotate  about  an 
axis  in  the  z-direction. 

A numerical  solution  for  the  systems  of  equations  governing  ■ 

each  problem  was  obtained  by  using  the  collocation  scheme  introduced  ' 

by  Erdogan,  Gupta,  and  Cook.  Results  were  computed  for  several 
values  of  the  geometric  parameter  y = c/h. 

The  accuracy  of  the  results  was  checked  in  several  ways,  and 
it  proved  quite  satisfactory  for  global  quantities  such  as:  1 

1)  Load  diffusion  from  the  block  into  the  elastic  half  space 
(vertical  displacement  problem) ; 

2)  Percentage  of  load  carried  by  the  sides  and  the  base  of  the 
block  (vertical  and  horizontal  displacement  problems) ; 

3)  Ratio  of  applied  moment  to  horizontal  load  (horizontal  dis- 
placement problem) ; 

4)  Moment  stiffness  (rotation  problem) ; 

5)  Normal  and  shear  stress  distribution  along  the  sides  and 
base  of  the  block,  excluding  the  vicinity  of  the  corners 
(+c,h)  (all  problems). 
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Since  the  corner  conditions  and  the  singular  behavior  near  the 
corners  (+c,h)  were  not  satisfied,  local  qualities  such  as  the  shear 
and  normal  stresses  are  only  valid  away  from  these  corners. 


x/c  .y 

|u  /u  j 
yb  y2 

y'  c/h 

0.3 

0.5 

1 

2 

4 

0.2 

0.3092 

0.0547 

0.0018 

0.0028 

0.0016 

0.4 

0.0154 

0.0119 

0.0101 

0.0011 

0.0010 

0.6 

0.2421 

0.0461 

0.0051 

0.0002 

0.0007 

0.8 

0.0258 

0.0132 

0.0092 

0.0047 

0.0015 

1 

0.0544 

0.0178 

0.0149 

0.0081 

0.0057 

TABLE  1 

Vertical 

Displacement: 

Comparison 

of  the 

y-displacements 

relative  to  the  midpoint  at  the  surface  y = 0 of  the  riyid 

block  (u  ) and  the  two- inclusion  (wl,  ) problems  for  several 
yb  y2 

values  of  y = c/h. 


c/h 


0.05 

-0.002 

-0.236 

1.232 

-0.061 

0.1 

-0.006 

-0.505 

1.550 

0.147 

0.2 

-0.020 

-1.026 

2.115 

0.372 

0.  3 

-0.036 

-1.568 

2.7  39 

0.664 

0.5 

0.002 

-2.091 

3.350 

1.356 

1 

1.123 

-0.083 

1.738 

2.913 

2 

6.424 

2.680 

-0.255 

6.656 

4 

27.685 

4.369 

-0. 959 

18.175 

"ABLE  2 

Rotation:  Variation  of  the 

to  the  rotational  stiffness 

contributions  of  the 

of  the  embedded  block 

stresses 

with  y. 
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Vertical  displacement  of  corner  (c,0)  relative  to  (0,0) 
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Fig.  1 Geometry  and  coordinate  system  for  a partially  embedded 
finite  rigid  block  undergoing  vertical  translation. 
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Fig.  2 Superposition  solution. 


Fig.  3 Geometry  and  coordinate  system  for  an  embedded  block 
undergoing  horizontal  translation  without  rotation. 
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y/h 


Fig.  5 Vertical  Displacement:  Comparison  of  load  diffusion  curves 
for  the  rigid  block  and  the  two- inclusion  problems. 
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Fig.  13  Rotation:  Normal  stress  distribution  at  y = h,  0<x<c 
for  various  values  of  c/h. 
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Fig.  15  Rotation:  Normal  stresses  acting  on  the  side  of  the 
block  for  0 < y < h,  for  several  values  of  c/h. 
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c/h 


Fig.  17  Rotation:  Stiffness  for  moment  applied  to  rigid  block 
(curve  Is  right  ordinate,  upper  abscissa;  curve  II: 
left  ordinate,  lower  abscissa). 


APPENDIX  A 


The  functions  K^,  K^,  L^i , and  L , i = 1,2,..,,  7,  are  given  as 


v - 2(h~  y)[(h -y)2  - (x-  s)2]  (K2-l)(h+y) 

Kni  ’ ’ [ (h-  y)2  + (x  - s)2 ]2  (h  + y) 2 + (x-  s)2 


(A . 1 ) 


+ 2<(h  - y)  [(h+  y)2  - (x-  s)2]  _ 8hy(h  +y)  [ (h  + y)2  - 3 (x  - s)2] 
[(h+y)2  + (x-  s)2]2  [(h+y)2  + (x  — s ) 2 ] 3 


(A.  2) 


K (x  ys)  - (x-  3)  f2U<  -2)(h-y)2  + k(x-  s)2]  (<2  + 1) 

Sr  ,y,S;  U [(h-y)2  + (x-  s)2]2 + (h+y)2  + (x  - s)2 

4k  (h+y)  2 8hy  [ 3 (h+y)  2 - (x-  s)2Ji 

[ (h+  y)2  + (x-  s)2]2  [(h+y)2  + (x-s)2]3  ^ 

„ , \ , Nf2[(K  + 2)(h-y)2  + <(x-s)2]  (k2+1) 

■ <»-  ,)V  + Th'i  y)n  d 

+ Mh  + Hi + 8hyl3(h  + y)2-  (x-  s)^]  i 

[ (h+  y)2  + (x-  s)2]2  [ (h+  y)2  + (x-s)2]^  ‘ 
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APPENDIX  B 

The  calculated  values  of  A,  B,  C,  and  D at  the  corner  (c,h) , 
defined  in  Chapter  VI,  are  given  in  this  Appendix  for  various 
geometries . 

The  values  given  in  Tables  B1-B3  are  based  on  the  numerical 
results  obtained  from  the  solution  of  the  problems  formulated  in 
Chapters  III-V,  respectively.  These  values  were  calculated  by 
means  of  a quadratic  extrapolation  based  on  the  three  points 
nearest  the  corner  (c,h). 

The  results  presented  in  Tables  B4-B6  are  those  obtained  from 
the  solution  of  the  considered  problems  when  three  of  the  corner 
conditions  are  incorporated  in  the  system  of  equations.  Therefore, 
these  results  satisfy  all  corner  conditions  exactly. 

As  stated  in  Chapter  VI  the  degree  of  accuracy  of  these  re- 
sults is  unknown.  However,  the  following  general  remarks  can  be 
made  in  association  with  these  results.  The  corner  conditions 
were  derived  by  an  asymptotic  expansion  of  all  singular  terms 
appearing  in  the  governing  equations!  this  includes  both  the  Cauchy- 
type  kernels,  and  the  terms  of  the  generalized  kernels  with  singu- 
lar contributions  near  the  end-points.  Moreover,  inclusion  of  the 
corner  conditions  in  the  system  of  equations  introduces  the 
proper  interaction  between  the  stress  discontinuities  in  the 
vicinity  of  the  corner  (c,h).  Also,  as  it  was  shown  previously, 
inclusion  of  the  corner  conditions  has  little  effect  on  the  global 


results . 
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Therefore,  it  seems  that  including  the  corner  conditions  in 
the  system  of  equations  besides  yielding  correct  global  results, 
will  also  improve  the  accuracy  of  the  results  in  the  vicinity  of 
the  corners,  provided  that  the  difficulties  discussed  in 
Chapter  VI  have  been  successfully  dealt  with. 


c/h 

A (1) 

B (1) 

C(l) 

D(l) 

0.05 

0.158 

-0.157 

0.203 

0.153 

0.1 

0.163 

-0.211 

0.229 

0.162 

0.2 

0.388 

-0.314 

0.278 

-0.024 

0.3 

0.599 

-0.316 

0.192 

-0.150 

0.5 

0.607 

-0.253 

0.089 

-0.168 

1 

0.551 

-0.295 

0.144 

-0.076 

2 

0.609 

-0.265 

0.079 

0.012 

4 

0.678 

-0.211 

-0.010 

0.039 

TABLE  B1  Vertical  Displacement:  Calculated  values  of 

A,  B,  C,  and  D at  the  corner  (c,h)  for  different 
geometries. 
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c/h 

A (1) 

B (1) 

C(l) 

D (1) 

0.05 

-0.107 

0.364 

-0.850 

0.146 

0.1 

-0.222 

0.569 

-1.026 

0.180 

0.2 

-0.356 

0.779 

-0.815 

0.134 

0.3 

-0.273 

0.661 

-0.435 

0.024 

0.5 

-0.136 

0.183 

0.199 

-0.097 

1 

-0.150 

0.011 

0.368 

-0.101 

2 

-0.240 

0.440 

0.029 

-0.103 

4 

-0.305 

0.604 

0.040 

-0.134 

TABLE  B2  Horizontal  Displacement:  Calculated  values  of 

A,  B,  C,  and  D at  the  corner  (c,h)  for  different 
geometries . 
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c/h 

A (1) 

B (1) 

C(l) 

D (1) 

0.05 

-3.630 

-29.868 

0.995 

-0.244 

0.1 

5.164 

-20.283 

1.470 

-0.349 

0.2 

7.262 

-16.335 

2.296 

-0.617 

0.3 

5.207 

-13.074 

3.127 

-0.866 

0.5 

2.303 

-7.217 

3.994 

-1.088 

1 

0.216 

-1.016 

2.083 

-0.536 

2 

-0.037 

0.037 

0.044 

-0.200 

4 

-0.022 

0.025 

-0.022 

-0.337 

TABLE  B3  Rotation:  Calculated  values  of  A,  B,  C, 
and  D at  the  corner  (c,h)  for  different 
geometries . 
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c/h 

A(l) 

B (1) 

C(l) 

D (1) 

0.05 

0.123 

-0.087 

-0.445 

0.316 

0.1 

0.144 

-0.102 

-0.346 

0.246 

0.2 

0.139 

-0.099 

-0.222 

0.157 

0.3 

0.134 

-0.095 

-0.167 

0.119 

0.5 

0.131 

-0.093 

-0.121 

0.086 

1 

0.158 

-0.112 

-0.097 

0.069 

2 

0.186 

-0.132 

-0.075 

0.053 

4 

0.182 

-0.129 

-0.049 

0.035 

TABLE  B4 

Vertical 

Displacement: 

Calculated 

values  of 

A,  B,  C,  and  D at  the  corner  (c,h)  for  different 
geometries  (using  corner  conditions) . 
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c/h 


A(l) 


B (1 ) 


C(l) 


D(l) 


0.05 

-0.004 

0.003 

0.015 

-0.010 

0.1 

-0.019 

0.014 

0.046 

-0.033 

0.2 

-0.081 

0.058 

0.129 

-0.092 

0.3 

-0.125 

0.089 

0.156 

-0.111 

0.5 

-C 

0.107 

0.139 

-0.098 

1 

1 

C 

1 

a 

K? 

0.129 

0.111 

-0.079 

2 

-0. 327 

0.232 

0.132 

-0.094 

4 

-0.569 

0.404 

0.152 

-0.108 

i 


TABLE  B5  Horizontal  Displacement:  Calculated  values  of 

A,  B,  C,  and  D at  the  corner  (c,h)  for  different 
geometries  (using  corner  conditions) . 
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c/h 


A (1 ) 


B (1) 


C ( 1 ) 


D (1) 


0.05 

-58.682 

41.616 

1.066 

-0.756 

0.1 

-28.960 

20.538 

1.393 

-0.988 

0.2 

-10.788 

7.651 

1.374 

-0.974 

0.3 

-5.418 

3.842 

1.220 

-0.865 

0.5 

-2.027 

1.438 

0.935 

-0.663 

1 

-0.316 

0.224 

0.385 

-0.273 

2 

-0.032 

0.023 

0.103 

-0.073 

4 

-0.015 

0.011 

0.130 

-0.092 

TABLE 

B6  Rotation: 

Calculated 

values  of  A, 

I,  C, 

and  D at  the  corner  (c,h)  for  different 
geometries  (using  corner  conditions) . 
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